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Motivated by the dark energy issue, a minisuperspace approach to the stability for mod- 
ified gravitational models in a four dimensional cosmological setting is investigated. Specif- 
ically after revisiting the f(R) case, R being the Ricci curvature, we present a stability 
condition around a de Sitter solution valid for modified gravitational models of general- 
ized type F(R,G,Q), G and Q being the Gauss-Bonnet and quadratic Riemann invariants 
respectively. A generalization to higher order invariants is mentioned. 

1 Introduction 

It is well known that recently it has been found strong evidence for an accelerated expansion of 
the universe, apparently due to the presence of an effective positive cosmological constant and 
associated with this acceleration there exists the so called dark energy issue (see for example 



The modified gravity models are pure gravitational alternative for dark energy (for a recent 
review and alternative approaches see [21 [3]). The main idea underlying these approaches to 
dark energy puzzle is quite simple and consists in adding to the gravitational Einstein-Hilbert 
action other gravitational terms which may dominate the cosmological evolution during the very 
early or the very late universe epochs, but in such a way that General Relativity remains valid 
at intermediate epochs and also at non cosmological scale. 

In the present paper we shall consider a large class of modified gravitational cosmological 
models defined in a Friedmann-Robertson- Walker (FRW) space-time and we shall focus our 
attention on the stability of the de Sitter solution. We recall that the stability issue is relevant 
in many contexts. For example, in the ACDM model it ensures that no future singularity will be 
present in the solution and within cosmological models, the stability or instability around a de 
Sitter solution is of some interest at early or later times. However, we remind that the inclusion 
of a cosmological term has to confront with the well known cosmological constant problem, an 
unsolved issue so far. On the other hand, as anticipated before, the modified gravity models 
may offer a quite natural geometrical approach in the spirit of the original Einstein idea. 

First we shall briefly revisit the class of models based on the action (here k 2 = 8ttGn, Gn 
being the Newton's constant) 
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where f{R) depends only on the scalar curvature R. We will consider pure gravitational theories, 
since we are interested in the dark energy sector, namely the property of the de Sitter critical 
point. The inclusion of ordinary matter can be done and it will be not treat here, even though 
its inclusion is important in reconstructing the expansion history of the Universe and probing 
the phenomenological relevance of the models (see for example the recent paper [U[5], where 
the case f(R) has been discussed). 

In a cosmological setting, the f(R) models have been introduced in [6j [7] and investigated 
in many papers (see [3]). Also the stability of the solutions has been discussed in several 
places |SllSlITDliniII2lII3lII3lHlIiaiISlirZl. To this aim, different techniques have been 
employed, including manifestly covariant and field theoretical approaches, where the gauge issue 
has been properly taken into account. All these investigations are in agreement with the following 
conditions which ensures the existence and the stability of the de Sitter solution: 

2/o = Rofa , (1-2) 

1<A- (1-3) 

where /', /" are the derivatives of f(R) with respect to R and /q, f{j are the derivatives evaluated 
at the value R = Rq. The first condition, Eq. (|1.2|) . determines the scalar curvature of the de 
Sitter solution, while the second one, Eq. (|1.3p . gives the condition for the stability around the 
de Sitter solution. In a more general quadratic theory, the stability has been investigated in 

PI EE]. 

The aim of this paper is to address the same investigation to a more general class of modified 
gravity, that is the one based on the action 

S = ^I dix ^F{R, G, Q) , (1.4) 

which depend on an arbitrary function of the scalar curvature R , Gauss-Bonnet invariant G 
and the quadratic Riemann invariant Q, the relation among them being 

G = R 2 -4P + Q, P = R^ V R?\ Q = R^paR^ ■ (1.5) 

The string-inspired scalar-Gauss-Bonnet gravity case F(R, G) has been suggested in Ref. [20] as 
a model for gravitational dark energy while some time ago it has been applied to the possible 
solution of the initial singularity problem [21j . 

The investigation of different regimes of cosmic acceleration in such string-inspired gravity 
models has been carried out in Refs. (2012211231121112511231231231123130!. In particular, in [28] 
a first attempt to the study of the stability of such kind of models has been carried out using 
an approach based on quantum field theory. 

The method we will use in order to study the stability for F(R, G, Q) models is a classical 
approach, which we call minisuperspace approach, and it is based on a Lagrangian formalism 
[311 l32"j I33j . inspired by the seminal paper [31], where, for the first time, quantum effects were 
considered. 

With regard to this, it is well known that one-loop and two-loop quantum effects induce 
higher derivative gravitational terms in the effective gravitational Lagrangian and early studies 
on instability for quadratic terms have been investigated in [35J. A particular case has been 
recently studied in [36] . 
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For general quadratic models, namely with linear dependence on R 2 , P and Q, the stability 
has been reported in p^8l f!9] . 

It should be stressed that the stability studied here is the one with respect to homogeneous 
perturbations. For the F(R) case, the stability criterion for homogeneous perturbations coincides 
with the inhomogeneous ones [9]. 

The content of the paper is the following. In Section II, we revisit the f(R) models and 
rederive the stability condition (|1.3p and then the same analysis is extended to the Gauss- 
Bonnet models in Section III and to the general case in Section IV. The paper ends with the 
conclusions. 



2 Minisuperspace approach for f(R) 



As already mentioned above, we shall deal with FRW isotropic and homogeneous models with 
spatial flat metric, namely 



ds z 



-N' z (t)dt z + a z {t)d 2 x 



(2.1) 



where t is the cosmic time, a(t) the cosmological factor and N(t) an arbitrary lapse function, 
which describes the gauge freedom associated with the reparametrization invariance of the min- 
isuperspace gravitational model. For the above metric, the scalar curvature reads 



R = 6 



iN 2 



+ 



a 2 



aN 



i 2 N 2 aN 3 



(2.2) 



As usual the 'dot' over the symbol means derivative with respect to time, so d = ^ and so on. 

If one plugs this expression in the Eq. (II. ip . one obtains a higher derivative Lagrangian 
theory. In order to work with a standard (first derivatives) Lagrangian system, we make use of 
a Lagrangian multiplier A and we write [3"2~] 



S 



1 

2^2 



dx / dtNa" 



f(R) ~ A 



R- 6 



+ 



a 2 



dN 



aN 2 a 2 N 2 aN 3 



Making the variation with respect to R, one gets 
df(R) 



A 



dR 



(2.3) 



(2.4) 



Thus, substituting this value and making an integration by part one arrives at the Lagrangian, 
which will be our starting point 



Lia, a, R, R, N) = -6 - 6 ^« + Na> 



N dR 



N d 2 R 



f(R)-R 



df(R) 



dR 



(2.5) 



Now, a and R are Lagrangian variables and iV appears as an "einbein" Lagrangian multiplier, 
reflecting the parametrization invariance of the action. Furthermore, the three equations of 
motion related to these three variables are not independent as we will show. In the analysis of 
the system we can use only two equations of motion and we can fix N(t) by the gauge choice 
(this corresponds to the choice of the cosmological time). 
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The first Equation of motion ^ = reads 

f'H + f'H 2 + Rf') = 0, (2.6) 

where for convenience we have introduced 

a an 

and we have chosen the gauge N(t) = 1. 

The conserved quantity related to the parametrization invariance is the energy computed 
with the standard Legendre transformation and it is vanishing on shell due to the equation of 
motion (|2.6f) for the einbein N. This is the energy constraint. 

The other equations of motion associated with the variation of R and a are respectively 

H + 2H 2 -- = 0, (2.8) 
6 

2f + if'H-2H 2 f'-^ + f = 0. (2.9) 

As we already said above, this latter equation is redundant. In fact, taking the derivative with 
respect to t of Eq. (|2.6|) and making use of Eq. (|2.8|) a direct calculation leads to Eq. (|2.9p . while 
equation (12. 8p is equivalent to (12.20 . 

The two equations (|2.6|) and (|2.8|) form a very simple autonomous system in the two variables 
R and H, namely 

H=--2H 2 . (2.11) 

6 V ; 

The analysis of stability is standard and consists first in finding the critical points Ro, Hq, which 
are defined by imposing R = and H = and then in investigating the stability of the related 
linear system around such critical points. In our case we have the solutions 

R = 12H 2 , (2.12) 

2/o-^o/o = 0, (2.13) 

which correspond to a de Sitter critical point with scalar curvature determined by the above 
condition. 

The linearized system around the de Sitter critical point reads 



Uh)= M { S h)' M =U -4flb ' ' <2 - 14) 



where M is the stability matrix. It is easy to show that the stability is ensured by the two 
conditions TrM < and detM > 0. The first one is trivially satisfied, while the second one, 
making use of Eq. (|2.13|) leads to the stability condition (|1.3p . 
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3 Minisuperspace approach for the Gauss-Bonnet model 



In this Section we shall generalize the previous approach to the modified gravitational model 
defined by a Lagrangian density of the type F(R,G). Here G is the Gauss-Bonnet invariant, 
which for the FWR metric reads 

24d 2 



G 



i 3 iV 5 



(aN - diV) 



(3.1) 



while R is given by Eq. (|2.2I) . 

In order to put the Lagrangian in a standard form, now we have to use two Lagrangian 
multipliers, so we start from 



A 



^2 J ' dx J dtNa 3 | F(R, G) 

(aN - aN 



R-e 



+ 



aN 



aN 2 a 2 N 2 aN 3 



24a 2 
a 3 N 5 



Making the variation with respect to R and G one gets 



A 



dF(R,G) 



F'r, 



dF(R, G) 



dR ' n 8G 

and by making an integration by part one arrives at the Lagrangian 

8d 3 dF' n 



(3.2) 



(3.3) 



L(a,a,R,R,G,G,N) 



6d 2 aF' R 
N 



G 



6aa 2 dF' R 



N 3 dt N dt 
+Na 3 (F — RF' R — GFq) 



(3.4) 



Of course, the couple of equations of motion related to the variations of R and G are equivalent 
to Eqs. ([231) and ((331) . which lead to 



R = 6(H + 2H 2 ), G = AH 2 {R-m 2 ) , 

while the equation ^ = gives 

24H 3 F G + 6H 2 F' R + 6HF G + (F - RF R - GF' G ) 



0. 



(3.5) 



(3.6) 



where we have chosen N(t) = 1 again and we have put H(t) = a(t)/a(t) as above. The conserved 
quantity is always the energy computed with the standard Legendre transformation and it is 
vanishing on shell. 

The equation of motion related to the variation of a is 

8H 2 F G + 2F' R + AHF' R + 16HF G (H + H 2 ) 



+ 



+ 6H 2 ) +F 



RF'r 



GF' G 



0. 



(3.7) 



Again, this latter equation is a consequence of the others. In fact, taking the derivative with 
respect to t of Eq. (|3.6|) and making use of Eqs. (|3.5j) . one obtains Eq. (|3.7|) . As a consequence, 
in order to deal with a first order differential autonomous system, we may use Eqs. (I3.5P and 
(|3.6|) . obtaining in this way 



R 



B(R,H) 
A(R,H) 



(3.8) 
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H = ^-2H 2 , (3.9) 
o 

where 

A{R, H) = QHF' R + 48# 3 i^ G + 96H 5 F GG , (3.10) 

B(R, H) = -F + (R- QH 2 )(F' R + AH 2 F G ) - 8H 2 (R - 12H 2 ) 2 (F'^ G + AH 2 F GG ) . (3.11) 

The critical points are defined by R = and H = 0. As a result 

Ro = 12flf , G = 24fl£ (3.12) 
and from B{Rq,Hq) = it follows 

F - G F G (R , G ) - Mk^Go) = Q (g ^ 

This corresponds to a de Sitter critical point with Gauss-Bonnet invariant determined by the 
condition (13.131) (see ref. [28J). 
Since 

19Bo -h rsui 
A- mT- Ho > (3 - 14) 

the linearized system around de Sitter critical point reads 

where M is the stability matrix and ^4o>-So are quantities in (|3.10p and (|3.11|) evaluated at Rq 
and Hq. Again, requiring TrM < ( again trivially satisfied) and detM > one obtains the 
stability condition 

i < — F Fr(RqiGq) ^ (316) 



Rn 



-^Rfl(-^O) Co) + ^RqF gr (R , Gq) + ^RqF gg (R , Go) 



It is easy to show that this condition, when F(R,G) = f(R) reproduces the relation f)l .3j) 
discussed in Section II. 

It has to be noted that the same result can be obtained starting from Eq. (|3.7p and reducing 
it to a first order autonomous system by the method discussed in [3], but in this case, the 
resulting stability matrix is a 3 x 3 matrix. 

We conclude the Section with an example. We recall that it is interesting to investigate 
Gauss-Bonnet models defined by 

F(R,G)=R + f(G), (3.17) 

since they may be relevant from a phenomenological point of view [291 [28]. in these special cases 
the conditions (|3.13|) and (|3.16|) read 

Gofo-fo = SH 2 . (3.18) 

(319) 
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As a particular example let us choose 

f(G) = aCP , (3.20) 
where a is a dimensional constant and 7 is a dimensionless parameter. Eq. (|3.18p gives 

2a( 7 - 1)G^ 1/2 = V6 (3.21) 



and this implies that 0(7 — 1) > 0, while Eq. (|3.19l) leads to 1 < 1/(27). Then the model is 
stable (around de Sitter) if the arbitrary parameters a and (3 satisfy both the conditions 



0(7 - 1) > 



2- 



> 1 



a < 
0< 7 <^ 



(3.22) 



All the other possible choices for the parameters a, (3 give rise to an unstable de Sitter solution. 
In particular, the model with negative parameters a and 7 is unstable. 



4 Minisuperspace approach for the F(R, G, Q) case 

Now we shall generalize the approach described in the previous section to a Lagrangian density 
of the type F(R,G,Q). The use of G instead of P simplifies the derivation. Here Q is the 
quadratic Riemann invariant, which for the FWR metric (|2.ip reads 



Q 



12 



i 2 N & 



(aN - aN 



2 12d 4 
+ 



x 4 7V 4 ' 



(4.1) 



In order to put the Lagrangian in a standard form, now we have to use three Lagrangian 
multipliers, so we start from 



^2 J dx J dtNa 3 I F(R, G,Q) - A 

Q - 



a - (aN - aN 



a 3 N 5 



R-6 



12 



a~ 



aN 2 



a?N 2 
2 



i 2 iV 6 



(aN - aN 



+ 



aN \ 

Up) 

12d 4 \ 
^iV 4 ) 



Making the variation with respect to R, G and Q, one gets 



A 



dF 



dF 



F'g, 



dF 



P 



F' n 



OR ~ n ' ^ dG ~ r dQ 

and by making an integration by part, one arrives at the Lagrangian 



(4.2) 



(4.3) 



L 



6a 2 aF' R 8a 3 dF' r 6aa 2 dF' R „ r „.„ „ . . 
K a K + Na 3 (F - RF' R - GF G - QF' Q ) 



N 



N 3 dt 



d 4 



N dt 

2 aa 2 N 2 



AT5 



2aaaN\ 
~N*~) ' 



(4.4) 



The equations of motion related to the variations of R, G and Q are equivalent to Eqs. 12.2 
(|3.ip and (|4.ip respectively and read 



R = 6(H + 2H 2 



G = 4H 2 (R-6H 2 ), 



Q = 12H A + -(R-6H 2 ) 2 . 
3 



(4.5) 
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Here we have put H(t) = a(t)/a(t) and we have choose the "gauge" N(t) = 1. The equation of 
motion for the variable N gives (in the gauge N(t) = 1) 

24H 3 F G + 6H 2 F' R + WF' R + (F - RF' R - GF' G - QF' Q ) 

+2AH(H + H 2 )F' Q + 12F' Q (2HH + 6HH 2 - H 2 ) = 0, (4.6) 

while the equation of motion related to the variation of a again is a consequence of the other 
ones and for this reason we do not write it. 

Now making use Eqs. 14. 5p and (|4.6p . we arrive at the following first order differential au- 
tonomous system again composed only by two equations, 

H = ^-2H 2 , (4.8) 
6 

where 

A(R,H) = 4HF f Q + 6H[F% R + 8H 3 (F% G + F% Q ) 

+ 1ZH\F'S G + 2F GQ + F% Q j\ + A x , (4.9) 

B(R,H) = -F + RF R + GF G + QF' Q -%H 2 F R + AH 2 F' Q (R-12H 2 ) + B 1 . (4.10) 

where A\ and B\ are trivially vanishing when evaluated at the critical points R = and H = 0. 
They are given by 

A 1 = 8H(R - 12H 2 ) (f%q + AH 2 F G q + |f£ q ) , (4.11) 

B x = -8H(R - 12H 2 ) 2 (f'^ g + AH 2 F GG - F% R ) . (4.12) 
As a result, one has 

Ro = 12H 2 , G = Q = 24i/ 4 (4.13) 
and from B(Rq,Hq) = 0, it follows 

Fo-%F G + Q F>)-^ = 0. (4.14) 



This generalizes the condition (|3.13p . 

The linearized system around de Sitter critical point reads 

(2)-"(S). "-{t^t*)' 

Aq = A(Rq,H\) being the quantity in (|4.9p evaluated at the critical point. Again the trace 
of the 2x2 matrix M is equal to the —3Hq < 0, while from detM > 0, we get the stability 
condition 

F' R + IRqFq 



1 < 



F'^ R + lF' Q + lR (F GR + F'^ Q ) + \Rl{F GG + 2F GQ + Fqq] 



(4.16) 
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This is the main result of our paper. If F does not depend on Q, the above condition reduces 
Eq. (I3,16p derived in Section III. 

Let us consider an example, namely the model defined by [37] 

F(R,G,Q) =R-fi 2 f(R,G,Q), (4.17) 

/(E,G)= K*TW' fo= ;/(i *°' Go ' g ° )= -rJ^TJTW (4 ' 18) 

where /i 2 ,a,/5, 7 are arbitrary dimensional and dimensionless parameters and n is an arbitrary 
integer number. Now from Eqs. (|3.13p and (|3.16p we respectively obtain 

^-(n + lV 2 /o = => / >0 => ( a + /3 + 7 r>0, (4.19) 

t 15 „ \ < , Unia,P,i) =2n(a + p + -f) + 2(3 + a. (4.20) 

w n (a,/3,7) 

Eq. (|4.19|) ensures the existence of the de Sitter solution while (|4.2(jp is necessary for the stability 
of such a solution. 

We immediately see that (|4,20p is not satisfied if all the parameters a, (3 and 7 have the 
same sign, but it is not satisfied also in some other cases. For example, it is never satisfied if 
both a + P + 7 and to n (a, P, 7) have the same sign. The cases in which a + P + 7 and u n (a, P, 7) 
have opposite signs require a more detailed analysis. 

As an example, we may take a = —7, 7 > 0, P > thus a + /3 + 7 = /3>0. Then, as soon 
as P < 2 (n+i) ' one nas a D e Sitter stable solution. 

We conclude this Section observing that the method we are dealing with is very suitable for 
generalization to the case where higher curvature invariants of thecurvature are present. An 
interesting example is the third order invariant Q3, given by 

Q3 = Rfj,uaf3R al3p<7 Rpa 

N 3 a 3 N 2 a 2 d Nha 2 S 3 d 6 . . 

a 3 N 9 a 3 N s a 3 N 7 a 3 N & a 6 iV 6 v ; 

We recall that this is the invariant which governs the two-loop ultraviolet divergences in pure 
gravity [351155]. 

If now we consider the special case F(R, G, Q, Q3) = f(R, G, Q) — bQs, then a straightforward 
calculation leads to the following relations for the existence and the stability of the de Sitter 
solution: 



^n-f {F' G + F' Q )-?^-^m-^ (1-22, 



1 < 



F' R + IRqF' + \rI 



R F RR + l F Q + l R o( F GR + F 'rq) + h R l( F GG + 2F GQ + F Qo) 



1 b r>2 
+ 3 K 



(4.23) 
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As an explicit example of this type, we consider the function 

F = R - 2A + a x R 2 + a 2 P + a 3 Q-bQ 3 , (4.24) 

where P is given by Eq. (jl.5p and in order to use Eq. (|4.22l it has to be replaced in terms of 
R, G and Q. From (l4~2"2"J) . we get 

J R -4A -^ J R|] = 0, (4.25) 

while condition (|4.23p leads to 

12 - bRl , . 

>0. (4.26) 



(3d + a 2 + a 3 )R + 

For A small, one has real positive solutions. When 6 = (the absence of the cubic term), the 
stability condition is 

a 2 + a 3 + 3ai > , (4.27) 



which is the result obtained in [19j . In particular, when a 2 = and 03 = 0, one obtains the well 
known result that models of the type R — 2Aq + a\R 2 are stable as soon as a\ > 0. 

As a final remark, one should note that for 6 > 0, besides the flat solution one can also have 
a de Sitter solution induced by quantum effects [51] starting with Ao = in the Lagrangian, 
the scalar curvature being Rq = 36/6. As it is well known, and evident from Eq. (I4.25p . this is 
not true for the quadratic case, which requires a non vanishing cosmological constant term Ao 
in the initial Lagrangian. In this case, the stability condition is satisfied when 

(a 2 + a 3 + 3 ai ) < -3Vb. (4.28) 

5 Conclusions 

In this paper we have presented the Lagrangian minisuperspace approach to the stability issue 
around a de Sitter critical point for a class of modified gravitational models depending on the 
Ricci scalar, Gauss-Bonnet and Quadratic Riemann invariants. It should be stressed that the 
stability studied here is the one with respect to homogeneous perturbations. For the F(R) case, 
the stability criterion for homogeneous perturbations coincides with the inhomogeneous ones [9]. 

With respect to other approaches, the method here presented has the advantage that the 
autonomous system one is dealing with, as far as the stability of de Sitter solution is concerned, 
consists always in two first order differential equations, simplifying considerably the stability 
analysis for arbitrary dependence on the chosen curvature invariants. 

The method has been first applied to the class of models depending only on the Ricci scalar R 
and the well known stability condition for this case has been recovered. Then the same approach 
has been applied to the scalar-Gauss-Bonnet models and its generalizations and a new general 
condition for the stability around a de Sitter solution has been found. 

With regard to this last issue, we would like to note that the conditions (I4.22|) (existence of 
the de Sitter solution) and and (|4.23p (stability condition around the de Sitter solution) are very 
general and they can be applied to a large class of Lagrangians, as shown in the last example of 
Sec, IV, where a third order invariant in the Riemann tensor has been investigated. 
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As far as the comparison with other methods is concerned, in reference |36j, a criterion of 
stability for a class of generalized modified gravitational models around de Sitter solution and 
described by a function of the kind W(R, Q — 4P) has been derived with a field theoretical 
approach. Since G = R 2 — AP + Q, our stability condition can easily be related to the one 
obtained within a theoretical approach. 

Furthermore, we also have to mention Ref. [30], where the model of Ref. [37] have been 
discussed. The model considered in that paper is described by the choice 

m,p,Q) = R- (aR2 + f P+cQr , (5.1) 

with a, b, c positive parameters and n a positive integer. The instability of the model is due 
to the appearance of a spin-2 ghost. As a result, an instability of the vacuum arises. In our 
notation, this model corresponds to the first model (|4.17|) we have studied in Section IV and we 
have 

a = a + ^, P = ~\ ^ = c+ \- ( 5 ' 2 ) 

When a, b, c and n are positive, our condition (|4,16p is not satisfied and, as a consequence, the 
model is unstable, in agreement with the conclusions reported in Ref. [40J. Thus, the instability 
is a generic feature of this class of modified models. Within a theoretical approach, this has 
been already noted in Ref. [4T] . 

Finally, we observe that the inclusion of the matter for a generic F might be non trivial, 
but, within our approach, can be done, for example, along the line of Refs. [H [5]. 
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